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T— I ■ Double parton correlations, having effects on the double parton scattering processes oc- 

o ■ 

' curring in high-energy hadron-hadron collisions, for example at the LHC, are studied in the 

^ I valence quark region, within constituent quark models. In this framework, two particle cor- 

, relations are present without any additional prescription, at variance with what happens, 

■ for example, in independent particle models, such as the MIT bag model in its simplest 
I— i' version. From the present analysis, conclusions similar to the ones obtained recently in a 

^ ^, modified version of the bag model can be drawn: correlations in the longitudinal momenta 

^ |. of the active quarks are found to be sizable, while those in transverse momentum are much 

CD ■ 

■ smaller. However, the used framework allows to understand clearly the dynamical origin of 
the correlations. In particular, it is shown that the small size of the correlations in trans- 

^ ^ verse momentum is a model dependent result, which would not occur if models with sizable 

, quark orbital angular momentum were used to describe the proton. Our analysis permits 

. therefore to clarify the dynamical origin of the double parton correlations and to establish 

which, among the features of the results, are model independent. The possibility to test 



O 

cn 



X 
J3 



experimentally the studied effects is discussed. 

I. INTRODUCTION 

Long time ago, the relevance of multiple hard partonic collisions taking place within a single 

hadronic scattering, the so called multiple parton interactions (MPI), has been addressed and 

I 

studied [1]. Although MPI are higher twist, i.e. they are suppressed by a power of Aqqy)/Q'^ with 
respect to single parton interactions, being Q the partonic center-of-mass energy in the collision, 
experimental evidence of these kind of processes has been obtained, already some years ago [^. 
MPI are of great relevance for LHC Physics, where they represent a background for the search of 
new Physics. It is not surprising therefore that a strong debate around MPI has been arising in 
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recent years, when several dedicated Workshops have been organized, starting from that illustrated 

n 

in Ref. [3]. In this work we concentrate on double parton scattering (DPS) processes. It is now 
understood that DPS contributes, e.g., to same-sign WW and same-sign dilepton production J-7| 
and represents a background for Higgs studies in the channel pp — )• WH — )• ivbb jsMll | . Evidence 



of the occurrence of DPS at LHC has been already established 



121 ]. Two comprehensive papers 



reporting the status of the art theoretical knowledge of these processes have recently appeared 
13, [l^. New signatures of DPS producing double Drell-Yan processes have been studied recently 
15l |. The dynamical origin of double parton correlations, having potential effects in DPS, in 

semi-inclusive deep inelastic scattering and in hard exclusive processes, has been the subject of 

recent comprehensive studies llGl l. 

n 

In the orig inalwork fll, the DPS cross section was written in terms of double parton distribution 
functions (dPDF) Fij{xi,X2, z±), describing the joint probability of having two partons with flavors 
hj = QyQydy longitudinal momentum fractions xi,X2 and transverse separation z± inside a hadron: 

d<T = ^ jd'^z±Fij{xi,X2,z±,n)Fkiix3,Xi,z±,n) 

i,j,k,l 

X aikixiX3^/s, ^J-)ajl{x2X4^/s, fj.) . (1) 

The partonic cross sections a represent the hard, short-distance processes, while 5 is a symmetry 
factor arising if identical particles are present in the final state. In Eq. ([1]), contributions due to 
diparton flavor, spin and color correlations are neglected, as well as parton-exchange interference 
contributions Q, Q, 17, 18]. These correlations are present in QCD and could be phenomenolog- 
ically relevant. 

Two main assumptions are usually made in DPS analyses: 

i) the dependence on the transverse separation and on the momentum fractions or parton flavors 
are not correlated: 

Fij{xi,X2,z^,fi) = Fij(xi,X2,/i)r(z_j_,/x) ; (2) 

ii) a factorized form is chosen also for the xi,X2 dependence: 

Fij{xi,X2,fJ-) 

= qi{xi,fj,)qj{x2,fJ.)0{l - xi- X2){1 -x\- X2Y , 

where q is the usual parton distribution function (PDF). 9(1 — xi — X2)(l — xi — X2)" introduces 
the kinematic constraint xi + X2 < 1, and n > is a parameter to be flxed phenomenologically. 
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dPDFs, describing soft Physics, are nonperturbative functions. Positivity bounds on dPDFs have 
been recently derived {22]. 

dPDFs cannot be evaluated in QCD but, as it happens for the usual PDFs, they can be at least 
estimated at a low scale, Qo ~ Aqcdi the so called hadronic scale, using quark models. 

The evolution of dPDFs, namely the way they change from a given momentum scale Qo to a 
given momentum scale Q > Qo is being a matter of debate jisl , Q, 3-21]. The resolution of 



this debate will be instrumental to relate not only data from different experiments but also model 
calculations at the hadronic scale to LHC data. In this way, the analysis of data involving DPS 
can be guided. 

A first important calculation of this kind is the one recently presented in Ref. 23], where a 



bag model framework for the proton [24] has been considered to evaluate dPDFs in the valence 
region. As in any model where the valence quarks carry all the momentum, in the bag model the 
initial scale /i has to be taken quite low. If the bag were assumed to be rigid, as in important 

the quarks would be 



early calculations, corresponding to the so-called cavity approximation 

independent and none of the relevant correlations described by dPDFs would be found. In Ref. 
therefore, a prescription to recover momentum conservation, leading to quark correlations in the 
bag, already applied in model calculations of parton distributions, has been used 
main outcome of Ref. 



20423]. The 

23l ] is that, in the modified MIT bag model, the approximation Eq. ([2]) holds 



reasonably well, but Eq. ([3]) is strongly violated. Problems with Eq. ([3]) had already been pointed 



out in Ref. 



3-: 



291-131]. 



The analysis of Ref. 



231 ] is retaken here in a constituent quark model (CQM) framework. CQM 



calculations of parton distributions have been proven to be able to predict the gross features of 



PDFs 



3d, 



331 ] ■ generalized parton distributions (GPDs) [34] and transverse momentum dependent 



parton distributions (TMDs 
respect to the approach of Ref 



Similar expectations motivated the present analysis. With 



231 ]. the non relativistic (NR) dynamics used here includes from 
the very beginning correlations into the scheme, so that, for example, the dynamical origin of 
the breaking of the approximation Eq. ([2]) is clarified. As it will be seen, this analysis helps to 
understand which, among the conclusions of Ref. 3], are model dependent. It is important to 



231 ] ■ is valid at the low 



stress from the very beginning that this kind of analysis, as the one in Ref. 
hadronic scale corresponding to the model and only in the valence quark region. For this study to 
be directly used in the LHC data analysis, having for the moment high statistics only for very small 
values of x, far therefore from the valence region, and corresponding to high momentum transfer, 
additional studies, in particular the introduction of the QCD evolution of the model results, are 
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necessary. This investigation, like the one in Ref. [23|, is therefore a first exploratory estimate, 
thought to understand the dynamical origin of the different momentum correlations and to guide 
experimental measurements accordingly. 

The paper is structured as follows. In the next section, the formalism to evaluate the dPDFs 
in a CQM is clarified and the main equations presented. In the third one, results are illustrated 
and discussed. Eventually, conclusions are drawn in the last section. 

II. CONSTITUENT QUARK MODEL DESCRIPTION OF DOUBLE PARTON 

DISTRIBUTION FUNCTIONS 

As already stated in the Introduction, we are going to use a CQM framework and therefore 
only correlations for two valence quarks of flavor qi, q2 can be studied. Let us work in momentum 
space, in terms of the Fourier-transformed distribution Fq-^^q^^xi, X2, k±), defined as follows: 

Fqiq2{xi,X2,k±) = f dz± e'"^^'^^ Fg^g^{xi, X2, z±) . (4) 



The formalism introduced in Ref. will be used. Color-correlated and interference dPDFs will 
not be considered, being Sudakov suppressed at high energies Q, [sgI- 

In a NR quark model, the color-summed dPDF for two quarks qi, q2 with polarization si, S2 in 
the proton, is: 



Fq^q2{xi,X2,k±) = J dkidk2 J dk[dk2nq^q^{ki,si;k2,S2;k[,si;k2,S2) 

where light cone components, given by = po i: for a generic 4-vector a'^ have been used 
for the quark and proton momenta and the two quark, spin dependent, off diagonal momentum 
distribution 

nq,q,{ki,si;k2,s2;k[,si;k'2,s2) = 3^ / d^g ^* (^i, si; ^2, ^2; ^3, ^a) A(1)A(2) A(1)A(2) 

S3 '' 

X ^ {k'i,si;k'2,S2;k3,S3) (6) 



has been introduced, using momentum space wave functions. In Eq. ([6]), the flavor and spin 
projectors, given by 

Puid){i) = 5 (7) 
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and 

- — 2 — ' 

respectively, for the particles 1, 2 appear. The factor 3 in Eq. ^ represents the number of valence 
quark pairs in the proton. It arises from the fact that any pair is equivalent to each other, since the 
proton ground state wave function is symmetric, once the color degrees of freedom are factorized 
out. 

By taking into account properly the corresponding flavor and spin projectors, any kind of spin 
and flavor correlations can be easily evaluated. For the sake of definiteness, in this paper we will 
concentrate on spin averaged dPDFs, which we call, for brevity, distributions qiq2- In a CQM 
framework, using translational invariance, they have therefore, according to Eqs. ([5]) and ([6]), the 
following form: 

qiq2{xi,X2,k_i) = dkidk2dkz'4>* i^i + ^,k2 -^,k^ Pq^{l)Pq.^{2) 

X ^\^ki-^M + ^M\6{ki + k2 + h) 5 (xi - ^ j (5 Us - ^ 



3 j dhdhr (^h + \^k2- ^ 4(1)^^92(2) 



X ^l,\^^-^,k2 + ^j5[x^-^y \X2 - -^J , (9) 

where the intrinsic wave function ip{ki, ^2) has been introduced and kj_ = \k±\. A similar expression 
was presented already in Ref. [14]; obtained as the NR limit of the general definition of the dPDFs, 
given in terms of light-cone quantized operators and states. This equation allows one to evaluate 
qiq2 in any CQM, by using the corresponding intrinsic momentum space wave functions. 
In our presentation, we will use two different CQM: 

i) a simple harmonic oscillator (HO) potential model, assuming an SU(6) spin-flavor structure, 
with the proton representing the ground state. In this case, the intrinsic wave function can be 
written: 

\N^S,/^)s = ijHki,k2) = ^ x SU{6) (10) 

where the spectroscopic notation \'^^'^^Xj)t, with t = A, M, S being the symmetry type, has been 
used and SU{6) is the usual SU(6) spin-flavor proton wave function (x refers to spin, $ to flavor): 

SU{6) = -^ixMS^MS + Xma^ma) , (11) 
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whose explicit expression is given in textbooks and not repeated here. The parameter o? 



muj 



of the HO potential is fixed to the value 1.35 fm~^, in order to reproduce the slope of the proton 
charge form factor at zero momentum transfer. Despite of its simplicity, it is perfectly clear that in 
this model dynamical correlations between the quarks are present from the very beginning, due to 
the ki ■ k2 term in the wave function, preventing it from being a simple product of two HO ground 
state wave functions. It is remarkable that this does not happen in the MIT bag_^model, whose 
simplest realization describes just independent, uncorrelated particles. In Ref. [23], correlations 
are introduced by properly modifying the simplest bag model scenario. In a potential model, we 
have them even in the simplest case. We think therefore it is very appropriate to study double 
parton correlations in this framework. 

ii) the CQM of Isgur and Karl (IK) [s^]. In this model, the proton wave function is obtained 
by adding a one gluon exchange contribution to a confining HO potential, along the line addressed 



already in Ref. 



; including contributions up to the 2fkj shell, neglecting a very small D wave 



contribution, the proton |A^) is given by the following admixture of states 



\N) = ■4){ki, k2,si, 82,83) = ag\N''^Si/2)s + a§'\N'^S[/2)s + «M|^^^'S'i/2)Af 



(12) 



The coefficients were determined by spectroscopic properties to be ag = 0.931, a§i = —0.274, 
391 ]. In this case, the parameter = muj of the HO potential is fixed to the value 



Am 



-0.233. 



1.23 fm 5L,in order to reproduce the slope of the proton charge form factor at zero momentum 
transfer 



391]. If ag = 1 and ag' = = the simple HO model is recovered. 
The formal ex pre ssions of the wave functions appearing in Eq. (jl2p . yielding the IK model, can 
be found in Ref. [3^, [4^, given in terms of the following sets of conjugated intrinsic coordinates 



R 



1 ^ 1 ^ ^ ^ 

-^(fi + f2 + fs) o K = -j={ki + k2 + ks) 



P 



1 
71 



(fi - r2) O kp 



1 



(ki - k2) 



A = -^(fi + f2 - 2f3) ^ kx = + ^2 



2fe 



(13) 



There are many good reasons to use the IK model to estimate DPCs. First of all, the IK is 



a typical CQM, successful in reproducing the low energy properties of the nuc. 
spectrum and the electromagnetic form factors at small momentum transfer j37l . 



eon, such as the 



39|. In studies of 



DIS phenomenology, IK based model calculations were able to describe the gross features of PDFs 



33l |. generalized parton distributions (GPDs) [34] and transverse momentum dependent parton 



distributions (TMDs) [35|]. 
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FIG. 1: The dPDF uu{xi,X2,ki_) in the HO 
model, Eq. p3|) . for X2 = 0.4 at four values of 



FIG. 2: The same as in Fig. 1, but in the IK 
model, Eq. (fT7| . 



From now on, for the sake of definiteness, we will consider qi = q2 = u, i.e., correlations between 
the two quark u in the proton. In this case, the normahzation of the dPDF, defined according to 
Eqs. dSD - ([m, is found to be 



/ 



dxidx2uu{xi, X2, k± = 0) = 1 , 



(14) 



corresponding, as it is natural, to the number of uu pairs in the proton. 

By inserting the wave function Eq. (jlOh in Eq. ([9]), one gets, in the HO model: 

uu"^{xi,X2,k^) = Ce-^ j dkidk2e-f^^'^^^'''^5 (xi - ^\ 5 (x2 - ^ 



(15) 



where C = (V3/(7ra^))^ and /(^i, k2,a) = 2{kf + A;! + " k2)lo? . The delta function is worked 
out considering the nucleon with mass M at rest and the quark with mass m ~ M/3 on mass shell, 



with energy k^ = \Jm^ + k'^. In this way one has, for i= 1,2: 



Xi 



K 
p- 



m? + kf — ki 
M 



(16) 



From this assumption, a little support violation arises, e.g., a little tail of the longitudinal momen- 
tum distribution is found for x > 1. From Eq. (jlSp . it is clear that uu depends on k±, only in the 
exponent of the Gaussian function outside the integral. By inserting the wave function Eq. ()12p 
in Eq. ([9]), one gets instead, in the IK model: 



uu^^ {xi,X2, k±) = Og uug + Og, uu§' + agag/ uugg/ + uuyi, 



(17) 
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where uug = uu^^ and 



C 

-e ^ 



[k± ■ (fci - k2)/a^f - 3[2f{ki,k2, a) + k'i/a'] + 9} 



X 5 { xi 



p 



6 X2 



^2 
P~ 



(18) 



C ^1- 

uuss'{xi,X2,k±) = -^e'^ 



(19) 



uuM{xuX2,k^) = ^e-^ I dkidk2e-f^'^'^^''^H\ 



kl + kl + Aki ■k2-^ 



1 r 
2 



kj_ ■ {ki - k2) 



+ ^T^^{kl-klf - k^-{ki + k2) \ g{kiM,kA_) 



X 5 [ xi 



(20) 



where 



g{ki,k2,k_i 



^(1 mi 1 - mi I 0)Y{^^inp)Y{_^^{nx 

mi 

^(1 m2 1 - m2 I 0)Y,^^{np,)Y,_^^{n, 



m2 

3 

167r2 



[cos 6'p/ cos + sin 6'p/ sm 0;^ cos {4>x — 4>p')\ 



X [cos 6'p cos + sin 0p sin cos (0A — (/'p)] • 



(21) 



In Eq. ([2T]) . the sohd angles Vt\, ilp, J^p/ are the ones defined by the vectors kpi = {ki — k2+k±)/V2, 
kp = {ki — k2 — k±)/V2 and kx = V^{ki +k2)/V2, so that the angles Op(pi), (pp^p')^ ^A) 0A appearing 
in Eq. (I2ip can be written therefore through the components of ki, k2, k_i_ according to the 
transformation laws Eq. ()13p . 

Looking at Eqs. ()18ll2ip it is clear that, in the IK model, with respect to the HO case, be- 
sides the one in the Gaussian function outside the integral, an additional dependence on k±_ is 
found. In particular, in Eq. (j20p . a peculiar k± dependence is introduced by the presence of orbital 
angular momentum of the quarks. Anyway, since the correlations under scrutiny are for unpolar- 
ized protons, no dependence on the orientation of k± remains. One can easily realize this result 
mathematically by rewriting the above expressions in terms of scalar products. 



9 



1.5 



1 



0.5 

1 



-0.5 



- 


fcx 


= 0.0 GeV 






= 0.2 GeV 




/■;\ 


= 0.4 GeV 


- 


lf\ — 

i \ 


= 0.5 GeV 




/ \ 
/ 






/ 






1 ' 






0.2 0.4 0.6 0.8 

Xi 



1 



FIG. 3: The ratio Eq. (|22]) in the IK model, 
for four values of k±_ . 

III. RESULTS AND DISCUSSION 



The results of the calculation are collected in Figs. 1-6. We have decided to show basically the 
same quantities presented in Ref. [9], because they describe all the features we want to emphasize 
and moreover they allow an easy comparison of the results of the different models. 

The dPDF uu{xi,X2,h_[_) evaluated in the HO model, Eq. (jlSp . for X2 = 0.4 at four values of 
k I , is shown in Fig. 1. These results are qualitatively and quantitatively similar to those in Ref. 
231 1 once correlations are added to the simplest bag model scheme. As already stated several times, 
in the present scheme correlations are present from the very beginning. The k± dependence found 
here is Gaussian, as a consequence of the HO wave functions used. Also in the model of Ref. 23 1 



a similar Gaussian behavior was found. In Fig. 2 , the same is shown in the IK model, using Eq. 
(jl7p . Although the results are qualitatively similar and even quantitatively not far from the HO 
ones, relevant differences are found, as illustrated in the next two figures. As a matter of facts, in 
the HO calculation, the approximation Eq. ^ perfectly holds: the k± dependence is only in the 
exponent of the Gaussian function outside the integral in Eq. (jlSp . In the IK model, it is not the 
case, as it is seen in Figs. 3, where the following ratio is shown 

nn(xi, X2 = 0.4, kx_) 



ri(xi,fc_L) 



(22) 



uu{xi = 0.4, X2 = 0.4, fcj 

which should be k± independent if the approximation Eq. ([2]) were correct. Actually, a little 
violation, comparable in size to the one found in Ref. 



231] is obtained. In the present calculation 
it is perfectly clear that the size of the violation of the approximation Eq. ([2]) is model dependent: 
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FIG. 4: The function uu{xi, X2,k± = 0) evalu- 
ated in the HO model for three different values 

of X2- 



FIG. 5: The ratio Eq. ^ in the HO model, 
for three different values of X2- 



it holds exactly if the ground state of the HO is taken, it is violated if mixture of states, with given 
0AM, are added. The fact that it is violated weakly is, again, a model dependent feature. In a 
model with stronger SU(6) breaking with respect to the IK one, strong violations, related to the 
0AM content of the proton, would arise. 

The extent to which the other approximation, Eq. ([3]), is violated, is shown in Figs. 4 and 5, 
in the HO model. In particular, the ratio 

uu{xi,X2, k± = 0) 



r2{xi,X2) 



U{X2) 



(23) 



shown in Fig. 5, should not depend on the choice of X2, if the ap proximation Eq. ([3]) were valid. 
Results in Fig. 4 are qualitatively similar to the ones in Ref. 23|: for example the maxima move 
towards lower xi values by increasing X2- In Fig. 5, this trend is confirmed but the maxima become 
bigger when X2 gets bigger, a behavior which is opposite to the one found in the corresponding figure 
of Ref. 



23l |. This can be understood thinking to the different behavior of the single particle PDFs 
appearing in the denominator. While one particle parton densities in the bag are not vanishing at 
X = and present a sizable tail at x < 0, the same does not happen here, where a little support 
violation is found only at a; > 1. The strong violation of Eq. ([3]) is also reported in Fig. 6. This 
figure has been drawn thinking that the HO framework is indeed a strongly correlated one. By 



defining dPDFs uu and PDFs u introducing a parameter f3, as follows 



uuf3{xi,X2, k± = 0) 



(4-/3) 



3/2 



Xi 



p- 



6 X2 



p- 



(24) 
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FIG. 6: The ratio Eq. (|26p . for: a) /? = 0, i.e., in an uncorrelated scenario; h) (3 — 0.25; c) (3 — 0.5; d) 
/3 = 1, i.e., in the correlated HO framework. 



it is easily seen that the HO model, a very strongly correlated one due to the HO potential, is 
recovered for /? = 1. (3 = represents instead an uncorrelated model, like, e.g., the genuine bag 
in the cavity approximation. Intermediate values of /3 may be similar to the modified bag model 



shown in Ref. 



23( 1 . This behavior is found indeed, as demonstrated by Fig. 6, where the ratio 



rp{xi,X2) = 1— — 1— (26) 

is shown. The ratio we obtain is actually rather flat for xi ~ 2:2 — 1, while some structures are 



seen in the calculation of Ref. 



23|. 



In closing this section, an important caveat is in order. We reiterate that the present estimate 
of quark correlations, as the one of Ref. {2^, showing important violations of the often assumed 
approximations Eqs. ([2]) and ([3]), is expected to be reliable only in the valence quark region, i.e. for 
X > 0.1, and at the low energy scale described by the CQM. Actual data from LHC are dominated 
by the very low x region (x < 10"'^), at a very high scale of momentum transfer. Any comparison 
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of this kind of results with data is instructive only once perturbative QCD effects are implemented 
by using QCD evolution to the experimental scale. Certainly, due to the coupling of quarks and 
gluons, correlations will be generated even if exact factorization is assumed at the low energy scale 
of the model. 

IV. CONCLUSIONS 

In many phenomenological applications, relevant, e.g., for the LHC experimental program, 
parton correlations arc neglected in modelling double parton distribution functions. The latter 
quantities are therefore factorized in two terms, one depending on the longitudinal momentum 
fractions of the two quarks, xi and X2, and another on the transverse momentum separation, k±, 
e.g. a sort of x — factorization is assumed. In addition to that, also an xi — X2 factorization 
assumption is usually made, e.g., the distribution on the longitudinal momenta xi,X2 is taken 
to be uncorrelated. This study aims at establishing to what extent constituent quark models 
support the two types of factorization, comparing the results with those of an important previous 
analysis, performed in a properly modified version of the simplest bag model. In this way, model- 
dependent features can be distinguished from model-independent ones. An analysis of double quark 
correlations has been therefore performed in the valence quark region, for spin averaged double 
parton distribution functions. Two different constituent quark models have been used, i.e., a 
simple non relativistic SU(6) symmetric model and its generalization, the model of Isgur and Karl, 
corresponding to a mild SU(6) breaking due to the exchange of one gluon. For the aim of studying 
correlations and of understanding their dynamical origin, the framework used here, correlated from 
the very beginning through the interquark potential, is to be preferred with respect to models where 
correlations are added through prescriptions modifying independent particle models, such as the 
simplest version of MIT bag model. The conclusions of our work are as follows. In both constituent 
quark models, the X1—X2 factorization is strongly violated, confirming a similar conclusion obtained 
already in the bag model framework. This feature seems therefore a model-independent one. In 
the SU(6) symmetric model, where the proton is described by a pure symmetric S wave, exact 
X — k± factorization is obtained. This features is mildly spoiled in the Isgur and Karl model, as 
it happens in the bag model framework. With respect to the latter, the present approach allows 
to understand more clearly the dynamical origin of the breaking of the x — k± factorization, and 
its relation to the orbital angular momentum of the quarks. It is also found that the amount 
of violation of the x — k± factorization is a model dependent feature. Certainly this effect is a 
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small one in the scenario discussed here, where the symmetric S wave is dominating. Models 
with more sizable contributions from components of the wave function with higher values of the 
orbital angular momentum would certainly show a stronger violation. In closing, it is important 
to stress that, for the results of model calculations of double parton distributions, obtained either 
in a bag model or in a constituent quark one, to be phenomenologically relevant to the LHC 
Physics programme, the evaluation of their QCD evolution is crucial. This will be the subject of 
further studies, together with the evaluation of the same observables using relativistic, Light-Front 
approaches. In this framework, successfully applied in Hadronic Physics in general (see, e.g., 41[) 
and for the calculation of PDFs in particular (see, e.g., 3]), important contributions from high 
values of the relative orbital angular momentum naturally arise. 
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